Spatial correlations in polydisperse, frictionless two-dimensional packings 



O 

(N 



O 



X3 
O 

o 



> 

(N 

o 



X 



C.B. O'Donovan, M.E. Mobius 
School of Physics, Trimty College Dublin, Dublin 2, Ireland 
(Dated: January 20, 2013) 

We investigate next nearest neighbor correlations of the contact number in simulations of polydis- 
perse, frictionless packings in two dimensions. We find that discs with few contacting neighbors are 
predominantly in contact with discs that have many neighbors and vice versa at all packing fractions. 
This counter-intuitive result can be explained by drawing a direct analogy to the Aboav-Weaire law 
in cellular structures. We find an empirical one parameter relation similar to the Aboav-Weaire 
law that satisfies an exact sum rule constraint. Surprisingly, there are no correlations in the radii 
between neighboring particles, despite correlations between contact number and radius. 
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Disordered packings of particles are the quintessential 
model for amorphous materials such as granular pack- 
ings [H, emulsions [2], wet foams [1] and glass formers 
While the contact number distribution and its av- 
erage near the random close packing density [11-0, [B-S] 
have been extensively studied in these systems, little is 
known about spatial correlations in the contact network. 
Various models that have recently been put forward to 
predict the density distribution of contact numbers 
and forces [lOj in random close packings implicitly 
assume the absence of such correlations. 

We address this question through simulations of a two- 
dimensional model system with polydisperse, frictionless 
soft discs. At the random close packing density (j)c, the 
discs just touch and have an average contact number (z) 
close to 4 as required for mechanical stability P, [ll], [13] . 
Due to disorder the individual contact numbers z are 
distributed according to some distribution P{z) that is 
typically non-Gaussian and depends on the polydisper- 
sity of the disc size distribution fi!]. Here, we investigate 
whether spatial correlations in the contact network exist. 

In our simulations we find that discs with many con- 
tacts favor neighbors with fewer contacts and vice versa. 
These correlations persist for all packing densities. This 
result is a direct analogue to the well known Aboav- 
Weaire law in the field of cellular structures which states 
that cells with fewer nei ghb ors are surrounded by cells 
with many neig hbors [HHll. We will show that our re- 



sults are in excellent agreement with a modified Aboav- 
Weaire law. Since geometrical constraints in the packing 
dictate that smaller particles have fewer contacts on av- 
erage , one may expect similar correlations for the size 
distribution in the packing, namely that larger particles 
are surrounded by smaller ones. Surprisingly, we find 
that the size of the central particle is uncorrelated to the 
average size of the contacting neighboring particles. 

We simulate the disordered packings by using Durian's 
soft disc model [sl , as implemented by Langlois et al. [3| ■ 
The discs have a harmonic repulsion proportional to their 
overlap and experience viscous tangential drag. We use 
1500 polydisperse discs whose size is normally distributed 



around the mean (r) with a variance (T2 = (0.304(r))^. 
The polydispersity allows us to access a wide range of 
contact numbers which is important to measure next 
nearest neighbor correlations. The discs are randomly 
placed in a periodic box at low packing fraction and then 
allowed to relax while their radii are slowly increased. 
The simulation terminates when the total elastic energy 
due to overlaps reaches a steady state at a predefined 
packing fraction. For each packing density, up to 10 dif- 
ferent packings are created to increase the statistics of our 
correlation measures. Upon reaching equilibrium, discs 
with fewer than three contacts (rattlers) are removed for 
the analysis of the contact network but are accounted for 
in the packing fraction. Contacts are defined as overlaps 
between discs. 

We study packings which range from the random close 
packing density (/)c up to (/) = 1.35. Note that the overlap 
area between bubbles is counted twice in the calculation 
of (j) in line with previous simulations of packings f]\ . In 
our simulations we find (jjc — 0.845 and the correspond- 
ing average contact number (z) = 4.07, which is close to 
the isostatic prediction (z) = 4 ^ . As shown in the inset 
of Fig[l] (z) increases approximately as 4 -I- 3.29y/(j> — (f>c 
close to the isostatic point, which is consistent with pre- 
vious results [1, [3]- We study packings up to (j) — 1-35, 
where the average contact number reaches (z) = 6. 

Figure [lJa)-[Hb) shows the distribution of the relative 
contact number, P{z— (z)), for different packing fractions 
and their respective variance /i2 = {z'^) — (z)'^, where 
(z*) — z^P{z). The shape of the distribution, which 
is not Gaussian, is independent of the packing fraction 
as evidenced by the collapse of P(z — (z)) onto a master 
curve, while ^2 varies only slightly. 

Next, we address the main result of our work - the 
correlations in the contact network of the packings at dif- 
ferent densities. Given a disordered packing of friction- 
less discs with a certain global average contact number 
(z), are the local contact numbers of neighboring parti- 
cles correlated? Here, we define neighbors to be discs 
in contact, i.e. discs that overlap. In order to quantify 
nearest neighbor correlations, we measure m(z), which 
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The Aboav-Weaire relation is a solution of the sum rule: 
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FIG. 1: (a) The contact number distribution P(z—{z)) for dif- 
ferent packing fractions 0: (x) 0.845; (*) 0.88; (V) 0.92; (+) 
0.96; (A) 1.0; (0) 1.07; (□) 1.14; (Q) 1-2; (>) 1.35. The inset 
shows the variation of the average contact number (z) with 
packing fraction 0. The red line corresponds to the square 
root scahng fit: (z) — 4 + zo\/(f> — (j>c, where zq = 3.29. (b) 
The variance of P{z), fi2, versus (2). 



is the average contact number of the neighbors of a disc 
with contact number z. 

This analysis is analogous to the pioneering work of 
Aboav [16] on polycrystals, where the same analysis was 
performed for the coordination number of cells in a cel- 
lular structure. He found that many-sided cells are sur- 
rounded by few-sided cells on average and vice versa. It 
has been noted by Weaire [17], that m obeys an exact 
sum rule which is independent of dimensionality.: 



(1) 



where z is the coordination number of the cells. This 
sum rule is based on a counting argument and one can 
show that the sum on the left hand side of the equation 
counts cells with z neighbors z times. In the absence of 
correlations in the coordination number between neigh- 
boring cells, m{z) is simply a constant (= m). It follows 
from the sum rule that in this case m = {z)+^2/{z) [l3| . 
In the context of cellular structures, this is referred to as 
a topological gas, although its existence is disputed [l^ . 



m = (z) — a + {{z)a + H2)/z, 



(2) 



where a is an empirical parameter . For a topological 
gas, a = — /i2/(z), but in a natural cellular structures 
such as dry foam a 
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Therefore, many sided 
cells have few-sided neighbors and vice versa. One can 
interpret this anti-correlation as a partial screening of 
the topological charge (z — (z)) by its nearest neighbors 
whose combined charge is z{m — (z)). Most 2D cellular 
structures, such as polycrystals and dry foams, obey this 
relation well [14i-il&] , with the notable exception of ran- 
dom Voronoi tessellations that exhibit deviations from 
this rule 

The counting argument that leads to the sum rule (Eq. 
P])) was originally developed for cellular structures but 
holds equally well for neighbors in a contact network of a 
disordered packing. The main difference is that in 2D cel- 
lular structures with threefold vertices, (z) — 6 13], while 
frictionless packings in two dimensions have an average 
contact number (z) greater than or equal to 4 depending 
on the packing fraction [l|, Q . 

The results for m are shown in Figl2l^a), where we plot 
(to — (z))z — fi2 versus (z — (z)) for three different pack- 
ing fractions. For the Aboav-Weaire law (Eq. ([2])) to 
hold we expect the data to follow a line with slope —a. 
Clearly, the data does not follow the uncorrelated pre- 
diction a = — /i2/(z), instead we observe spatial corre- 
lations: discs with few contacts are surrounded by discs 
with many contacts and vice versa. Another key result 
are the deviations from purely linear behavior, especially 
at higher packing fractions. 

In order to account for this non-linearity, we expand 
the Aboav-Weaire law in terms of the moments of the 
contact number distribution such that it still satisfies the 
sum rule 



(m - (z))z - = - 2^ Ci{z^ 



(3) 



where the c^'s are arbitrary constants. If q = for i > 1, 
one recovers the usual Aboav-Weaire law with ci = a. In 
order to fit our data it proved sufficient to only make C2 
non-zero, which leads to 



= (z) -bz + 



fj.2{l + b) + bizY' 



(4) 



where b = 02- This is a one parameter fit, similar to the 
Aboav-Weaire law. However, mz is now quadratic in z, 
instead of linear. As shown in FiglH^b), Eq. @ captures 
the non-linearity well and leads to a much improved fit 
compared to Eq.([2])- Including higher order terms in the 
expansion (Eq. 0) does not improve the fit significantly. 
The inset of Figl2][b) shows the decrease of the parameter 
b with (z), which means that the screening of topological 
charge decreases as the isostatic point is approached. 
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FIG. 2: Nearest neighbor correlations of the contact num- 
ber, (a) (m — {z))z — ^2 versus {z — (z)) for three different 
densities: (x) 0.854; (A) 1.0; (O) 1.2. The error bars are 
standard deviations from the mean. The lines are fits to Eq. 
Q. The dotted lines correspond to the uncorrelated predic- 
tion a — —fi2/{z). (b) [(m — {z))z — /i2]/& versus (z^ — (z^)) 
for all densities (same symbols as in Fig[lla)). The dotted 
line corresponds to the slope —1. The inset shows the fit 
parameter b as a function of (z) . 



We would like to stress that the analogy between cor- 
relations in the coordination number in disordered fric- 
tionless packings and cellular structures is not obvious, 
since these systems are governed by different local and 
global constraints. Although polydisperse packings can 
be tesselated into a cellular structure f3|, not all faces of 
a cell correspond to contacts, therefore the existence of 
correlations in packings does not follow naturally from 
similar correlations in disordered cellular structures. 

There is also a correlation between contact number and 
particle size in polydisperse packings Larger parti- 
cles have more contacts on average, since one may fit 
more particles around them on average. Figure |3] shows 
the average contact number of a particle with radius 
r, {z\r) = where P{z\r) is the conditional 

probability of a particle of radius r to have a contact 
number z. It is well described by a linear relation 

{z\r)^{z){l + j{r-{r))). (5) 
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FIG. 3: The average contact number {z\r) for a given particle 
radius r at different packing fractions </!>: (x) 0.845; (V) 0.92; 
{□) 1.14. Lines are fits to Eq.([5]). The inset shows the fit 
parameter 7 as a function of 0. 



Since {z\r) is constrained by the equality 
{z\r)P{r)dr = (z), there is only one empirical 
fit parameter 7, which varies little with cj) (FiglH inset). 
However, this linear relationship does break down at low 
r, since {z\r) > 3. A similar result exists in the field of 
cellular structures and is known as Lewis' law [l3. [20|. 

Given the nearest neighbor correlations in the coordi- 
nation number (Eq. ^) and the correlation between size 
and contact number (Eq. (O), are smaller particles sur- 
rounded by larger ones, similar to observations in cellular 
structures [22J ? 

In order to study radii correlations, we measured 
{Rnnlf), which is the average radius of neighboring discs 
in contact with a disc of radius r. Figure S] shows 
{Rrm\f) / {^) versus r / {r) for (pc- No correlations are ap- 
parent and the result agrees well with the uncorrelated 
prediction i?„„ which is discussed below. 

The counting argument that leads to the Weaire sum 
rule for the coordination number can also be applied for 
the radii. Analogously, the average radius of the neigh- 
bors (i?n„|r) needs to satisfy the following relation 

1*00 poo 

/ {Rnn\r){z\r)P{r)dr ^ / r{z\r)P{r)dr. (6) 
Jo Jo 

The left hand side of the equation amounts to an integral 
over the disc radii r weighted by {z\r). In the absence of 
correlations {Rnn\f) is a constant (= Rnn), and we have 

_ r{z\r)P{r)dr _ r{z\r)P{r)dr 

J^nn — poo, I \ -r,/ \ 1 ^ / \ ' ^ ' 

Jo {z\r)P{r)dr {z) 

Substituting the empirical relation for {z\r) (Eq.lH)), we 
find that i?„„ = {r){l+^a2/ {r)'^)^ which is slightly larger 
than (r) and varies little with (inset FigHJa)). At 0c, 
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weaker but does not vanish and is well described by a 
modified Aboav-Weaire law. This result is a direct ana- 
logue the topological screening observed in cellular struc- 
tures. Nevertheless, the physical origin of the screening 
parameter b remains unclear much like the a parameter 
in the Aboav-Weaire law • 

We want to emphasize that the counting argument that 
leads to the sum rule for the contact number and radii are 
valid in any dimension. Therefore, one may expect sim- 
ilar correlations in three dimensional packings as well as 
in frictional packings 2l[ . It remains to be seen whether 
these correlations depend on the preparation history of 
the packing, which is known to have an influence on (jjc 
0. 
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FIG. 4: Nearest neighbor correlations of the radii (a) {R„n\r} 
versus r/{r} for (p = 0.845. The dotted line corresponds to 
the uncorrelated prediction i?„„ (Eq.((7}). The inset shows 
Rnn versus (z). (b) Rnn / {Rnn\r) versus r/(r) for all densities 
(same symbols as in Fig[lja)). 



we obtain i?„„ = 1.042(r), which is consistent with our 
results from Fig|4l[a). At higher packing fractions shown 
in Figlll^b), {Rnn\r) / Rnn remains constant and close to 
1. Only for high and low r, slight deviations due to low 
statistics are observed. 

While the absence of correlations in the radii for neigh- 
boring particles may be expected given our preparation 
procedure where particles are placed in the box at ran- 
dom, it is surprising in the light of the two correlations we 
have measured. Namely, the correlations between con- 
tact numbers of neighboring particles (Eq. (U)) and the 
correlation between size and contact number (Eq. (O). 
The reason for this counter-intuitive result is that the 
relationship between the average contact number m and 
the corresponding average radius Rnn does not follow the 
linear relation (Eq. ([5])) for a single disc [2^. 

Although we have only shown results for packings with 
normally distributed radii, similar correlations are ob- 
served for other polydispersities such as bidisperse dis- 
tributions ^2^. 

In conclusion, we studied polydisperse, frictionless 
packings at various packing densities. We find that discs 
with many contacts are surrounded by discs with few 
contacts and vice versa. As the isostatic point is ap- 
proached, the screening of topological charges becomes 
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